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We study straight vortices with global longitudinal currents in the Bogomolny limit of the Abelian 
Higgs model with two charged scalar fields. The model possesses global SU{2) and local electromag- 
netic L'^(l) symmetries spontaneously broken to global f/(l) group, and corresponds to a semilocal 
limit of the standard electroweak model. We show that the contribution of the global SU{2) cur- 
rent to the vortex energy is proportional to the total current squared. Locally, these vortices carry 
also longitudinal electromagnetic currents, while the total electromagnetic current flowing through 
a transverse section of the vortex is always zero. The vortices with high winding numbers have, in 
general, a nested pipelike structure. The magnetic field of the vortex is concentrated at a certain 
distance from the geometric center of the vortex, thus resembling a "pipe." This magnetic pipe is 
layered between two electrically charged pipes that carry longitudinal electric currents in opposite 
directions. 

PACS numbers: 03.75.Mn, 03.75.Lm, 12.15.-y 
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I. INTRODUCTION 

Certain extended topological defects are able to sup- 
port longitudinal currents. A well-known example of 
such a defect is the Abrikosov-Nielsen-Olesen vortex so- 
lution [l| in a C/(l) X U{1) model 0: the corresponding 
vortex string may carry both bosonic and fermionic su- 
perconducting flows of large magnitude. Loops of these 
superconducting strings may have various cosmological 
consequences [3]. 

The extended topological structures with longitudinal 
currents appear also at much smaller scales in the context 
of condensed matter physics. For example, vortices that 
carry longitudinal superfluid flow in a superfluid '^He-A 
were proposed theoretically in Ref. [J]. Later, certain 
signatures of these stringy objects, known as w vortices, 
were found in nuclear magnetic resonance experiments [5| 
(see Ref. [y| for a detailed review). 

The symmetry patterns and the zoo of the topologi- 
cal defects in the superfluid ^He have a lot in common 
with the corresponding properties of the standard elec- 
troweak model in particle physics j7|] . Thus it is not sur- 
prising that one of the most interesting realizations of the 
current-carrying strings in the field theory context ap- 
pears in the standard model of electroweak interactions. 
This model supports various (embedded) topological de- 
fects including well-known Z- and W vortices @, Q . In 
a semilocal limit of the model (characterized by a spe- 
cial value of the Weinberg angle, 9w = '"'/S) these vortex 
solutions are known to carry persistent longitudinal cur- 
rents of isocharge associated with the global symmetry 
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subgroup [To]. The word "semilocal" indicates that the 
symmetry group of the model is a product of the local 
electromagnetic group, C/e.m.(l), and the global isosym- 
metry group, SUi(2). In this limit the original non- 
Abelian gauge group decouples from the scalar fields, and 
the resulting theory is basically the Abelian Higgs model 
with two complex scalar fields. 

Beyond the semilocal limit (i.e., at any value of Wein- 
berg's angle) the current-carrying solutions were con- 
structed in Refs. [Ill, [l2l and evidence of their pertur- 
bative stability was found (la. We refer an interested 
reader to the reviews [ij] and [ij] for an extended discus- 
sion of the vortex solutions in the standard electroweak 
model. 

In the semilocal limit the vortices with nonzero longi- 
tudinal isocurrents exist only in type-II regime tlfl] be- 
cause the energy per unit vortex length of such vortices 
diverges as the system approaches the boundary between 
type-II and type-I superconductivity. This boundary is 
known as the Bogomolny limit [la|. General vortex so- 
lutions without longitudinal currents were described in 
the Bogomolny limit of the semilocal model in Ref. [ig] . 
The (iso) charged vortex solutions with persistent longi- 
tudinal currents were found in the same limit in Ref. [17{ . 
We continue this line of investigation by studying in de- 
tails the current-carrying vortex solutions in the Bogo- 
molny limit of the Abelian with two Higgs fields. We 
concentrate on the current-carrying vortices with multi- 
ple winding numbers because these vortices - contrary to 
the elementary vortices with a unit winding number [l7[ 
- do have a finite energy per unit vortex length. 

Our study is also motivated by the (theoretically an- 
ticipated and/or experimentally found) existence of var- 
ious systems with multiple-component order parameters 
such as the already mentioned superfluid '^He f6'|, two- 



component Bose-Einstein condensates [T^ , two-band su- 
perconductors that includes the weh-studied system of 
MgB2 [laLiliquid metahic hydrogen [2^ , two-component 
plasmas [2l| etc. 

The structure of this paper is as follows. In Sec. [TT] 
we describe the Abelian model with two charged scalar 
complex fields. This model admits vortex solutions to the 
classical equation of motion, which are drastically simpli- 
fied in the Bogomolny limit. Section HII] is devoted to the 
illustration of various multivortex solutions in this limit. 
We concentrate on the vortices with high winding num- 
bers which have a "pipelike" structure: the magnetic field 
in such vortices is nonzero only in a thin region located 
at a certain fixed distance from the geometrical center of 
the vortex. In Sec. IIVI we introduce a longitudinal wave 
along the vortex and check that this wave corresponds 
to a current-carrying solution that satisfies the original 
(second-order) classical equations of motion. Finally, in 
Sec. |V] we calculate the energy per unit length of these 
vortices. The last section is devoted to our conclusion. 



According to Eq. ^ both components of the scalar 
field carry the same electric charge e. 

The dimensionality of the symmetry group (|4]) indi- 
cates that there are four conserved Noether currents. The 
electric (super) current corresponds to the rotations in the 
t^e.m.(l) direction. 
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J^,^p'j^. = -[^W^^-c.cl 



2i 



(7) 



For the sake of further convenience we introduced here 
the reduced electric current: 



if, = Jf./p'^ 



(8) 



The three other Noether currents can be grouped into 
the isovcctor field 



Kf^ = 2eJf^n — p n x df_in 



<j)tf(j) 



n' = l. 



(9) 
(10) 



II. TWO-COMPONENT MODEL 

The Lagrangian of the two-component Ginzburg- 
Landau model is 
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1 



C = ^-F^^.F^^'' + - \D^^' - A (|$|2 - n^ 



(1) 



where F^,i, = d^A^ — d^A^ is the field strength tensor of 
the Ue.m. (1) gauge field A^ and D^ = 9^-|-ieA^ is the cor- 
responding covariant derivative, and |<I>p = p^ = $t<j) 
The scalar field $ has two complex- valued components. 



■!> = 



Pie" 
P2e'- 



s''^! cos 61 
^1 e'^^sinB 



(2) 



where we introduced the parametrization in terms of the 
two phases (^1.2 G [— 7r,7r), two condensates pi^2 with 



P = \Ip\ + pI . 



and/or the angular variable 6* € [0,7r). 
The model ([T]) possess the symmetry 

G = 5(7/(2) x;7e.„,(l). 



(3) 



(4) 



where SUi{2) group rotates the Higgs fields in the iso- 
topic space. 



SUj{2) : I * 



fj, ^ Afj, , 



(5) 



with U G SUi{2), and the Abelian C/e.m.(l) group cor- 
responds to the local electromagnetic symmetry which 
affects the phases of the scalar field ([2]), 



Ue.mXl) : 



A.. 



A^ - id,M^) . 



(6) 



The isocurrent K^ corresponds to the isospinor rotations 
SUi{2). The unit isovector ft - that is constructed with 
the Pauli matrices r = (r^, t'^ ,t^) - transforms in adjoint 
representation of the SUi{2) subgroup ([5]). The Abelian 
gauge transformations (|6]) leave both the isocurrent (|9]) 
and the vector PUI) intact. 

It is convenient to make the rescaling, 

$^77$, A^^T^A^, x^^^, (11) 

er] 

so that all variables appear to be dimensionless. Then 
the classical equations of motion of the model ^ are 

5^Jm-+p'j.-^9'^/m. =0, (12) 
-d^p + ^pd^nd'^n + pj^f - ^ap (p^ - l) = , (13) 



where we denoted a — A/e^, and 



(14) 
(15) 



are the strength tensors for the (reduced) electric current 
Jm' ^'i- ®' ^^'i f°^ th'5 vector field n, Eq. pOj) . respec- 
tively. Instead of the gauge- variant field A^ we use below 
the gauge-invariant electric current j^ as an independent 
variable. The change from gauge- variant variables A^^ 
and <& to the gauge- invariant fields j^, n and p, is com- 
monly used in this model [22, 1231 . 

Equations ([T^ and P^ are supplemented by the con- 
servation laws of the electric current ([7]) and the isocur- 
rent ^, respectively: 



du.J^, 



0, 



d^K^ 



0. 



(16) 



III. PIPELIKE VORTEX SOLUTIONS IN 
BOGOMOLNY LIMIT 

In this section we discuss known static straight vor- 
texhke solutions to the classical equations of motion (fT2|) 
and ([T3| . These vortices do not carry longitudinal cur- 
rents so that the corresponding solutions depend only on 
two-dimensional coordinates xi and X2- Since all vec- 
tors/tensors involving (at least one of) the components 
/i = 0, 3 are zero, the system becomes effectively two- 
dimensional. 

In a special limit, 



A 



(17) 



the classical equations of motion p^ and ([T5)) reduce to 
a simpler system of differential equations: 
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X dx 
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nxdjn]^0, (18) 
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0, (19) 



where x — \x\. The expressions for the electric current ([7]) 
and the isocurrent ([9]) are simplified as well: 



J,;- 



jpdjp, 



K, ^ ~e,jd,{p''n) . (20) 



The special choice of couplings p7|) is known as the Bo- 
gomolny limit 15]. The vortex solutions of the "semilo- 
cal" Abelian two-Higgs model ~ that were first found 
in Ref. M - were subsequently studied in this limit in 
Refs. [9', 'l6]. A general solution to the classical equa- 
tions of motion in the Bogomolny limit (|17p has a nonzero 
condensate ([3]) at its origin and in general is dubbed as 
"skyrmion" (for a detailed review see Ref. |13|'). 

In the Bogomolny limit the magnetic fiux of the vortex 
is quantized: 



T: 



where 



/27r 
d X B — — k , 



fceZ, (21) 



B = F,2=3i2-^h2^^{l~p'). (22) 

2e 2 

is the magnetic field of the vortex. The first integration 
in Eq. (PT|) is taken over the spatial contour Coo with an 
infinitely large radius. The last relation in Eq. (P^ is 
supported by equations of motion (|19p and Eq. ((20|) . 

The string tension (i.e., the vortex energy per its unit 
length) is quantized as well (we restore r) for a moment): 
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^l^l^vrry^l^l. 



(23) 



where i = 1,2 are two-dimensional indices in the trans- 
verse (with respect to the straight vortex) plane. 



The validity of the quantization (1211) can be checked 
in the angular parametrization of the scalar field ^. 
The energy functional (|23|) contains the scalar derivative 
squared. 



+ ^sin'20[a^(^2-<^i)] 



(24) 



For the energy to be finite, a vacuum state is to be re- 
alized at the spatial infinity. This means that all four 
terms in Eq. (P^ should be zero: the variables 6 and 
p 7^ should be independent of the coordinates, and the 
electric current, 

Jj = p[eAi + cos^ Odi^pi + sin^ 9diip2) , (25) 

should be vanishing at spatial infinity. Moreover, if at 
spatial infinity the angle 9 docs not take the specific val- 
ues, Ooo 7^ and 6oo 7^ 7r/2, then the two phases of the 
Higgs fields should be equal, 1^9 1 = Lp2- We then recover 
the quantization of the magnetic flux (j2ip . where the 
number k gets interpreted as common winding number 
of the phases (pi and ip2- If either ^^^ = or 6*00 — 7r/2, 
then If I ^ tp2 and at large distances the winding number 
is carried either by ipi or by (^2, respectively. 

Using a stereographic projection one can parametrize 
the vector n, Eq. (|10|) . by a complex function u ~ u{z) 
with z — xi + ix2'. 



ni + in2 



2u 
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l + \u\^ 



na 



1 



(26) 



Then the flrst equation of motion (|18|) is satisfled auto- 
matically provided u is a meromorphic function |16j . 

For a vortex located at the origin of the two-dimen- 
sional plane, the following boundary conditions are to be 
implied: 



lim n^ (0,0, 1)'' 



limn= (0,0,-1) 



(27) 



where the flrst relation is the choice of the vacuum state. 
In terms of the function u the conditions (|27|) read as 

u(^ = 0) = oo, m(z = oo) = 0. (28) 

The suitable meromorphic function is 



where k is the mentioned winding number and 
a = eryi?vort 



(29) 



(30) 



is a dimensionless parameter which defines the size i?vort 
of the vortex core. This parameter is equivalent to the 
parameter ^0 of Q and to the parameter |go| of p^ . 

Substituting the solution (|29|) into Eq. ([26]) and, sub- 
sequently, into Eq. p^ . one gets an expression for /12. 
Then, introducing a normalized coordinate 



e 



X 

a 



Rv 



(31) 



we rewrite the second equation of motion (1191) as follows: 



ld_ 



d 



a% 2 



1 



(^2fe + 1)2 



(32) 



This equation fixes the behavior of the scalar conden- 
sate p as a function of distance from the vortex center ^, 
Eq. ([ST]) : for fixed vorticity k, Eq. ([ST]) : and for the fixed 
size of the vortex core a, Eq. (150]) . 
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FIG. 1. Condensate p at the distance r from the vortex center 
for various sizes (thicknesses) of the vortex cores i?vort ~ a/er). 
The vortices have vorticity numbers fc = 1 (top), fc = 2 (mid- 
dle), and fc = 10 (bottom). 



carrying a single unit of the flux, fc = 1 (the plot at 
the top); a double-flux vortex, fc = 2 (the plot at the 
middle); and the high- vorticity solution, fc = 10 (the plot 
at the bottom). One can observe that while the single- 
flux solutions p — p{C) are always monotonic, the higher- 
flux solutions are not. A nonmonotonic density variation 
but for a different kind of vortices in a two-component 
model was also found in Ref. [2^. The solutions with a 
large-sized core (a 3> 1) and with a large flux (fc = 10 
in our example) have a well-recognized minimum around 




The solutions of Eq. ([32|) can be found numerically. In 
Fig. [1] we show the condensate p, Eq. (jH]), for a vortex 



vort 



FIG. 2. Condensates p, pi and p2 in high-vorticity (fc = 10) 
solutions with two sizes a = erjRvoTt of the vortex cores. 

The condensates for these vortex solutions are nonzero, 
p 7^ 0, both at their centers, r = 0, and at spatial infini- 
ties, p -> oo. As for the individual condensates, the pi 
condensate is vanishing at the center, while the p2 con- 
densate is zero at spatial infinity. This particular pat- 
tern of the condensates in the vortexlike solution is quite 
common for many-component systems. It appears in par- 
ticle plwsics in the context of the standard electroweak 
model [13], and in the condensed matter physics: the cor- 
responding structures were experimentally observed in a 
two-component condensate of ^7Rb atoms [i8|. We plot 
the individual condensates of two solutions with a = 4 
and a = 10 for the fc = 10 vortex in Fig. [21 

In general, the magnetic field (122]) of a vortex with a 
multiple vorticity is concentrated outside the center of 
the vortex. In Fig. [3] we show the density plots of the 
magnetic field in the transverse (a;i,X2) plane for fc = 10 
vortices of different sizes i?vort ■ One can clearly see that 
the maximum of the magnetic field (darker regions) is 
distributed in a (generally) cylindrical thin region of the 
radius R ~ i?vort = a/erj. Thus, the magnetic field of the 
vortex is concentrated, basically, within a thin pipe of 
the radius i?vort- The magnetic field inside and outside 
the pipe is zero. A possible existence of similar solutions 
with pipelike ( "helical" ) structure in this model was dis- 
cussed in [2^. The annular vortex solutions that some- 
what resemble our pipelike solutions were also discussed 
inRefs. dHH. 

The pipelike structure is well pronounced for large 
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FIG. 3. Strength of the magnetic field (|22p in high-vorticity 
{k = 10) solutions in the transverse plane. Various sizes of the 
vortex cores -Rvort = a/erj are shown. The spatial coordinates 
are given in units of erj. 



(a ^ 1) vortices with high-vorticity numbers. In a dif- 
ferent context, vortices with high winding numbers were 
also studied in the standard electroweak model in [23| ■ 

Let us now consider the spatially-transverse compo- 
nent of the conserved isovector current ^ which cor- 
responds to the flow of the global SUi{2) charge. The 
diagonal (in the isospace) component of this current is 



,(3) 



d 



Krm.a)^-^^- 



2k 



1 



2k 



p^{t,k,a) 



(33) 



In the polar coordinates of the transverse plane, the ra- 
dial component of the transverse isocurrent ([55)) is zero, 
Kr =0. The angular component of this current is plot- 
ted in Fig.|3]for a few vortices with the large (fixed) wind- 
ing number fc = 10. The angular component exhibits a 
familiar pipelike structure with two peaks at moderate 
values of the vortex size a. These peaks are merged into 
one peak at larger values of a. 



IV. LONGITUDINAL CURRENTS 

The vortices may host certain waves that propagate 
freely along the vortices. The simplest way to introduce 
such a wave is to generalize the meromorphic function 
u [given, for a static solution, by Eq. (125]) ] as follows 
(remember that z — xi -\- 1x2)'- 



t(xi,X2,X3,a;4) = (^)''e^("^°-P=^^^) 



(34) 



Here the wave frequency ui = po and the wave momentum 
P3 are two unknown parameters. According to dill) , in 
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FIG. 4. The spatially transverse diagonal isocurrent p3p for 
vortices of various sizes a with a large winding number k — 10. 



Eq. p4p the wave frequency w and the wave momentum 
P3 are given in units of erj. 

According to Eq. I^E^ the ansatz (IMl) corresponds (up 
to a gauge transformation) to the following behavior of 
the scalar doublet ([5]): 



$ 



(35) 



Here r and ip are the polar coordinates in the transverse 
plane and pi_2 are certain functions of r. The form of 
Eq. ([55)1 is essentially the same as in Refs. [13] and [l7| . 
The ansatz p4p for the meromorphic function u{z) was 
used in a different context in a 0(3) model in Ref. [28|. 
In analogy with the previous section, we substitute 
the ansatz (p4|) into the classical equations of motion, 
(IT^ and (IT51) using Eq. (I^5|) . The new equations also in- 
clude the electric density, Jq, and the longitudinal elec- 
tric current J3, that are treated in these equations as 
the unknown independent functions. However, one can 
immediately figure out that the second-order differential 
equations p^ and ([T5)) are self-consistent provided that 

a; = ±P3, Jq{x)^^Jz{x)= jp^ . (36) 

The first relation in Eq. ([5S)) indicates that the longi- 
tudinal current is carried by massless waves. The second 
relation means that the regions with nonzero electric cur- 
rent density are always electrically charged. The signs in 
Eq. (|36p distinguish between two possible directions of 
the wave propagation along the vortex. 

The presence of the longitudinal currents does not in- 
fluence the behavior of the condensate p = p{r) because 
Eq. ([5^ is valid in the Jo 7^ case. The transverse 
components, p — 1,2, of the electric current J^ are not 
modified by the presence of the longitudinal current as 
well, and Eq. (PH)) remains unaffected. As for the longitu- 
dinal components, the (reduced) electric charge density 
j = JqI (? , Eq. ([S|) - or, equivalently, the longitudinal 
current density, Eq. (I36p - is determined by the follow- 
ing differential relation: 
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4ujk^i 
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(^2fe + 1)3 



ir - !)■ (37) 



A. Longitudinal electric current 

The total longitudinal electric current 13°' = 
/dj;ida;2 J3 streaming through any transverse plane 
of the vortex, and the total electric charge p*°* = 
/ dxida;2 J3 per unit vortex length are both zero. Indeed, 
an integration of the left- and right-hand sides of Eq. ((37)) 
over the spatial coordinates of the two-dimensional trans- 
verse plane gives us the vanishing result: 



±h 



d^rj{r)p'^{r) = 0. 



(38) 



Here we used the second equation in (|36p and noticed 
that C^j(C) — ^ in the limits ^ ^ and ^ — > 00. 
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FIG. 5. Magnetic field B (the dashed lines) and the longitu- 
dinal electric current J (the solid lines) in the A; = 10 vortex 
as the function of the radius r. The red/thin (blue/thick) 
lines correspond to the vortex of the size a = 4 {a — 10). 

Despite the total electric current flowing through any 
transverse plane is zero, the local structure of the vortex 
in terms of electric currents is, however, nontrivial. In 
Fig. [5] we plot the density of the electric current and the 
strength of the magnetic field as functions of the distance 
to the vortex with the vorticity k — 10. We show these 
quantities for two vortices with the sizes a = A and a — 
10. In both cases the electric current changes the sign 
at a certain nonzero distance Rq from the center of the 
vortex, 



J{Ro) = 0, 



< Rq < 00. 



(39) 



For a typical vortex the alternation of the current takes 
place at the value of the radius, that is approximately 
equal to the size of the vortex core, Rq ~ i?vort- As 
one can see from Fig. [Sj it is the very radius where the 
magnetic field takes its maximum. The density of the 
current is positive in the inner core of the vortex, r < 
i?vort and negative in the outer space, r > Rvoit- 

In Fig. |6] we show a two-dimensional plot of the dis- 
tribution of the electric current in the transverse plane 
of the fc = 10 vortex of the size a = 10. The current is 
concentrated in the two nested narrow pipes that have 




FIG. 6. The density of the longitudinal current - and, equiv- 
alently, the electric charge density - in the transverse plane 
of the fc = 10 vortex of the size a = 10. The inner (red) 
pipe corresponds to the electric current flowing in the posi- 
tive direction, while outer (blue) pipe denotes the density of 
the negatively directed current. There is almost no current 
inside the inner pipe and outside the outer pipe and also in 
between these two pipes. 



relatively large diameters. The electric current is almost 
zero outside these pipes. Notice, that due to Eq. (pG]) the 
current-carrying pipes are also electrically charged. The 
electrically charged pipes enclose the pipe of the mag- 
netic field, which was visualized in the lower right plot 
of Fig. [3l The interior and exterior of the vortex are 
oppositely charged, as illustrated in Fig. [S) 

The distribution of the transverse electric currents in 
the transverse plane have a qualitative similarity with 
the distribution of the longitudinal currents. Fig. [3l The 
transverse currents - defined by the first relation in 
Eq. ([20]) - rotate in opposite directions inside and out- 
side the pipe of the magnetic field. The transverse cur- 
rent density vanishes approximately at the radius of the 
magnetic pipe, similarly to the density of the longitudinal 
currents. These properties are the natural result of the 
nonmonotonic behavior of the condensate p(r). Fig. [T] 
Physically, the outer current (say, circulating the vortex 
clockwise) generates the magnetic field along the z axis 
in an "outer" tube of a certain radius. The inner current 
(that rotates in the counterclockwise direction) leads to 
appearance of the magnetic field in the opposite direc- 
tion in the "inner" tube. Since the inner and outer tubes 
have different radii, the resulting profile of the magnetic 
field has a maximum at a certain nonzero radius, thus 
forming a pipelike structure. 

The total positive and negative longitudinal current, 
/(^) = ±/^, can be calculated by integration over, say, 
the interior of the vortex core: 



/^|/(±)| 



dVj(r), 



(40) 



r^Ro 



where we take, for definiteness, J{r) > 0. The radius 
Rq is defined as a nonzero finite distance from the vortex 
center at which the electric current vanishes (j39|) . The 



integrated current (l40l) is shown in Fig. [7] as a function 
of the vortex size a for a few fixed vorticities. For smah 
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FIG. 7. The net positive/negative current (140 |l carried by 
interior/exterior of the vortex with the winding number k vs 
the vortex size a — e'qRvort- The inset shows zoom in on the 
region with a < 1. 



vortices, a < 1, the total positive/negative current is very 
small. However, the inner and outer pipes of the thick 
vortices, a > 1 may carry rather strong electric current. 
Moreover, the larger the transverse size of the vortex the 
stronger is the electric current |29|. 

Summarizing this section, we visualize the internal 
structure of the pipelike vortex in Fig. |51 The magnetic 
field has a pipelike profile shown by the green surface 
which is sandwiched in between two other pipes. The 
inner and outer pipes carry the electric currents in pos- 
itive (shown by the red surface) and negative (the blue 
surface) directions, respectively. Since these electric cur- 
rents in these pipes have both longitudinal and transverse 
components, the three-dimensional structure of the elec- 
tric currents resembles spirals which are also visualized 
in Fig. E 




FIG. 8. A graphical representation of the structure of a 
thick (a 2> 1) vortex with a high winding number (|fc| S> 
1). The pipelike magnetic field - shown by the green surface 
in the middle - is layered in between positive (inner) and 
negative (outer) electrically charged pipes (shown by the red 
and blue surfaces, respectively). These electric pipes carry 
helical electric currents in opposite directions. The arrows 
mark directions of the helical electric currents. 



B. Longitudinal isocurrent 

Now let us now consider the longitudinal components 
of the conserved isovector current ©■ According to 
Eqs. (|26l) and (IMl) . the off-diagonal components of the 
isocurrent, K^ and Kjj, , are oscillating in time and 
space because they depend linearly on periodic functions 
cost? and sini?, where -d = uj{xo ± x^) — ip. The diag- 
onal components are time-independent quantities. The 
longitudinal isocurrent is 



iff (e;fc,a) = -2c^ 



fP{C, k, a) 



2k 



c 



2fc 



^2*^ + 1 



+(i-e''^)j(C;fc,a) 



(41) 



The isocharge density, Kq, is linked with the longitudinal 
isocurrent density, K^, in a manner of relation (1361) for 
the electric charge/current density: Ksix) — ±Ko{x). 
Both the isocharge density and the longitudinal isocur- 
rent are proportional to the frequency parameter uj. 

The longitudinal isocurrent (PT|) is shown in Fig. [HI 
This current has a two-peak structure for all studied val- 
ues of the core size a. The peaks are centered near the 
values of the radius that approximately correspond to the 
extremes of the longitudinal electric current, plotted in 
Fig. [S] for a = A and a = 10. 

In addition. Fig. |9] demonstrates that - unlike the or- 
dinary electric current - the isocurrent density is either 
positively or negatively valued. Thus, the net isocur- 
rent streaming along the vortex is nonvanishing and the 
vortex has a nonzero total isocharge. These global prop- 
erties are qualitatively similar to the characteristics of 
the solutions with low winding numbers that were found 
in Refs. [l^l and [13. 
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FIG. 9. The longitudinal isocurrent (|41|l for vortices of vari- 
ous sizes a with large winding number k = 10. 

The absolute value of the total longitudinal isocurrent 
which is carried by the vortex. 



^tot _ 



dViff^(r) 



(42) 



is shown in Fig. \TU\ as a function of the vortex core size 
a for a few values of the winding number k. The behav- 
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FIG. 10. The net longitudinal isocurrent (|42|l vs the vortex 
core size a = erjRvort for various winding numbers k. 



ior of the new isocurrent P^ is similar to the net pos- 
itive/negative electric current, Fig. [T] The total isocur- 
rent is very small (at fixed uj) for small values of the 
vortex core size a. At a particular vortex size a, the 
isocurrent starts to grow quickly as a function of a. Then 
it reaches a maximum, and, unlike the positive/negative 
fraction of the electric current, the isocurrent slowly di- 
minishes as the function of the vortex core size a. In 
general, the larger is the winding number fc, the smaller 
is the net isocurrent. 

Notice that the total longitudinal electric charge car- 
ried by the solutions vanishes exactly ([55)1 because this 
current is generated by the local U{1) symmetry ([6]). The 
net longitudinal isocurrent is not zero (j42p since it carries 
the charges corresponding to the global SU{2) subgroup 
of the model ^. 



V. ENERGY OF VORTICES 

A. Contribution of longitudinal current 

In this section we calculate the energy of the vortices 
that is carried by the longitudinal currents in the Bo- 
gomolny limit. The energy momentum tensor of the 
Abelian two-Higgs model (UJ is given by the following 
equation. 






(43) 



where g^i, is the metric tensor. In the spherical isovector 
representation (fTO|) the energy density reads as follows: 



Too = T^idopf + T^{d.P? + ^{d^nf + P-{dSf 



2 9 

P ,2 , P -2 



-Hp'-iy + ^j'o + ^j^ 






(44) 



A ^"-3 o /i. 



We integrate Eq. (|44|) over the transverse coordinates Xi 
and X2j imply the Bogomolny limit (J17p . and simplify 
the expression using both integration by parts and the 
corresponding equations of motion (IT51) and (|18l) . As a 
result, we get the vortex energy per unit vortex length: 



a{k,a) = d^xTQQ ~ TT\k\ + x{k,a)io'^ 



(45) 



The first term in the right-hand side of Eq. (|45|) corre- 
sponds to the energy of the vortex in the absence of the 
electric current (|23p . The second part is the contribution 
of the electric current to the vortex energy. This con- 
tribution is proportional to the wave frequency square, 
w^ or, equivalently, to the square of the total isocurrent 
carried by the vortex (|42p , or to the total electric current 
in, say, the inner charged pipe of the solution (|40l) . 

The dependence of the energy on the current is gov- 
erned by the dimensionless coefficient x that can be sep- 
arated into three parts: 

X{k,a) =xi(k,a) + X2{k,a) + xsik) ■ (46) 

The first term in Eq. (j46p is expressed via the conden- 
sate p. 



Xi{k,a)^2TTa' / d^ 



p+V'(^;fc,a) 







(^2fc + 1)2 



(47) 



while the second term can be determined via the electric 
current density j 



X2{k,a)=8TTk' /de 





1^ 



(^2fe + 1)3 



2fc-l 



-j{^;k,a) 



J=l 



(48) 



The third term can be calculated exactly: 



x3ik) = 87Te /de 



? 



4fc-l 



27r|fc| 



ie" + 1)' 



(49) 



The energy of the vortex (H5|) saturates the Bogomolny 
bound and it thus corresponds to the lowest possible en- 
ergy for fixed parameters fc, uj and a. Expressions ((45])- 
(gni) correspond to the lower bound on energy calculated 
previously in Ref. [U^ . 

The first integral (j47]) diverges for fc — 1. Indeed, 
at large transverse distances ^ ^ I the condensate 
is distance independent, p ~ 1, and the integral in 
Eq. (|47)) becomes logarithmically divergent at large dis- 
tances, J "" ^~\i^ ~ log^max- Thus, fiuite-sized vor- 
tices (a ^ 0) with a minimal winding number, fc = 1, 
cannot carry the longitudinal currents in the Bogomolny 
limit. However, the vortices with higher vorticities, 
fc ^ 2, can support the longitudinal currents in the this 
limit. 

In the Bogomolny limit the energy of the currentless 
vortex - given by the first term in Eq. (j45p - is indepen- 
dent of the vortex size a. Unexpectedly, in the presence 



of a longitudinal current the energy becomes dependent 
on the vortex size parameter a. In Fig. [TT]we plot x as a 
function of the vortex size a for the same set of vorticities 
that was already selected for Fig. [71 Notice that despite 
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FIG. 12. The energy coefficient x as a function of the (nor- 
malized) absolute value of the electric current flowing in, say, 
the inner charged pipe (|40p . 



FIG. 11. The coefficient x that determines the contribution 
of the electric current into the energy of the vortex (|45p . The 
inset shows a zoom in on a low-a region. 

the fact that the curves in Fig. [7] and in Fig. [TT] look sim- 
ilar on a qualitative level, they cannot be superimposed 
on each other by simple scale transformations. 

According to Fig. [TTl the excess of energy provided by 
the electric current is a monotonically increasing function 
of the vortex size at relatively large size a of the vortex 
core. A zoom in on the low-a region - shown in the inset 
of Fig. [TT]- reveals a very shallow minimum at the fixed 
size of the vortex 



0.8 , or i?v 



0.8 
er] 



(50) 



for all studied vorticities fc ^ 2. 

Thus, in the Bogomolny limit a vortex with a fixed 
winding number fc > 1 and a fixed frequency parameter 
u! y^ would tend to change its transverse size to the 
value (|50l) that corresponds to the global minimum of 
the vortex energy. According to the inset of Fig. [71 the 
electric current at a ~ amin is small but nonzero (we 
remind that the current can be made arbitrarily large by 
increasing the free parameter oj). 

The vortex core tends to stabilize its transverse size 
at the value (|5D|) for which the nested-pipe structure of 
the vortex interior cannot be resolved visually. This fact 
is seen in the behavior of the condensate p of the a = 1 
vortex with the vorticity A: = 10 (lowest panel of Fig. [TJ, 
or one can alternatively look at the magnetic field pro- 
file of the same vortex in Fig. [31 However, the minimum 
is very shallow so that a small perturbation ~ for exam- 
ple, induced by a shift of parameters towards type-II su- 
perconductivity - may shift the stability point to higher 
values of the vortex core radius a. 

The dependence of the vortex energy per unit vortex 
length on the value of the electric current in one of the 
pipes / is shown in Fig. [TH At very low values of the 
current / the energy has a very shallow minimum which 



is visible in the plot of the k = 2 solution. At higher 
values of the current the energy is a monotonically rising 
function. 



B. Rearrangement of vortex lattice 

As we have already discussed, in the Bogomolny limit 
the energy per unit vortex length of an elementary {k = 
1) vortex is infinite at any finite w, or, equivalently, at 
any nonzero value of the longitudinal isocurrent. In the 
absence of the external isocurrent the total energy of a 
many-vortex configuration does not depend neither on 
positions of vortices nor on their sizes in this limit. Now, 
in a thought experiment, let us apply an external isocur- 
rent 1^'°* along a configuration of an even number of 
parallel elementary vortices. Because of energy consider- 
ations it is clear that the elementary vortices would tend 
to merge together and form a set of, at least, double 
(fc = 2) vortices that have a finite energy at the nonzero 
isocurrent. Then each of the fc = 2 vortices would tend 
to adjust its transverse size a in such a way that the total 
energy is minimal at the fixed value of the isocurrent [30| . 
Thus, the presence of an external isocurrent should gen- 
erally lead to rearrangement of the vortex ensembles. 

A nonzero isocurrent (I4ip corresponds to a combina- 
tion of two electric currents - carried by up and down 
components of the order parameter - for which the total 
electric current is zero. Although a realization of such 
an isocurrent in an experimental setup is an open issue, 
our result suggests that the vortex lattice should be re- 
structured under the influence of the external isocurrent 
that is parallel to the external magnetic field. Also, our 
result shows that the multiple merging of the elementary 
vortices under the influence of the external isocurrent 
may lead to a formation of a vortex with a high winding 
number that exhibits the exotic "nested-pipe" structure. 

We expect that our results are rather generic and they 
should not be limited strictly to the particular limit P7)) 
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of the couplings. If the system is sufficiently close to 
the Bogomohiy limit, then we expect that an external 
isocurrent would induce a rearrangement of the vortices 
in this system. 

One of the well-studied examples of the two-band su- 
perconductor is MgB2 [19]. In external magnetic field 
the vortex lattice in MgB2 has quite specific geometrical 
patterns 31|, that are imprinted by elementary vortices. 
Such an inhomogeneous distribution of vortices in mul- 
ticomponent systems was first predicted theoretically in 
Ref. [321 ■ It would be interesting to check the effect of 
the external longitudinal isocurrent on the properties of 
the vortex lattice in this real material, which is rather far 
from the Bogomolny limit ([T7| . 



VI. CONCLUSION 

We found pipelike vortex structures in an Abelian 
gauge model with two scalar condensates. This model 
possesses global SUi{2) symmetry (isosymmetry) as well 
as local Ue.m. (1) symmetry corresponding to the Maxwell 
electromagnetism. The symmetry group is broken spon- 
taneously to a global U{1) subgroup. The model pos- 
sesses vortices that are linelike topological defects which 
carry magnetic fiux. We found that the vortices with 
high winding numbers and with large sizes of the vor- 
tex cores have the pipelike shape of the magnetic field: 
the magnetic field is concentrated at a certain distance 
from the geometric center of the vortex, thus resembling 
a pipe. 

For the sake of simplicity we were working in the Bogo- 
molny limit of the model in which the classical equations 
of motion are drastically simplified. 

We show that nonelementary vortices in the Bogo- 
molny limit of the model are able to support the lon- 
gitudinal isocurrents that carry the isocharge along the 
vortices. The isocharge corresponds to the conserved 
Noether charge with respect to the diagonal component 



of the global SUi{2) subgroup. The isocurrent is pro- 
portional to the difference between the electric currents 
carried by the upper and lower components of the order 
parameter. The current-carrying vortex also possesses a 
nonzero density of the isocharge localized in the vicin- 
ity of the vortex core. The elementary current-carrying 
vortices with unit vorticity have infinite energy per unit 
length, while the energies of the current-carrying vortices 
with multiple winding numbers are finite. 

The total vortex energy is given by a sum of two terms. 
The first contribution is the standard Bogomolny term 
that is proportional to the magnetic flux inside the vor- 
tex. The second term comes due to the presence of the 
longitudinal current. This term is proportional to the 
total isocurrent squared. The dependence of the vortex 
energy on the isocharge current may lead to potentially 
observable effects like, for example, rearrangement of the 
positions of vortices in the vortex lattice of a two-band 
superconductor. 

The vortices with large transverse sizes and with high 
winding numbers have a nested pipelike structure. In 
addition to the mentioned "magnetic pipe" that carries 
the magnetic flux of the vortex, the vortex also possesses 
two electrically charged pipes of larger and smaller radii. 
These "electric pipes" are oppositely charged and they 
also carry longitudinal electric currents of equal strength 
in opposite directions so that the net electromagnetic cur- 
rent flowing through a transverse section of the vortex is 
always zero. The magnetic pipe is always layered between 
the two electric pipes. 
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